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Abstract: We discuss the role of nucleon and its excited state poles in the twist expansion 



of the nucleon structure function moments. We find that the nucleon pole contribution was 
overestimated in previous analyses by a factor of two. Inclusion of this missing factor makes 



the duality appear for all moments at least down to Q = 1 GeV . For resonance poles the 
time reversal invariance together with unitarity demand at least four of them in the forward 
Compton scattering amplitude. These poles as well as the elastic pole can be propagated 
separately through the standard Operator Product Expansion derivation for DIS. This part 
of the amplitude gives a coherent, positive, higher twist contribution, which can be singled 
^ ■ out of the total and compared with the remaining one given mostly by threshold effects. 

A comparison of the estimated resonance contribution to the data on structure function 
moments allows to test large-Q 2 behavior of proton-resonance transition form-factors. 
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1. Introduction 

Investigations of the nucleon structure with electromagnetic probes at large momentum 
transfers allow for a simple treatment within the parton model. In fact, in this kinematics 
the inclusive lepton-nucleon scattering cross section can be related to the parton momen- 
tum distributions inside the nucleon. Latter do not follow from pQCD and they are subject 
to a measurement or Lattice simulations. Since structure functions themselves cannot be 
compared to pQCD predictions one can instead extract quantities calculable in the theory, 
moments. The Mellin transform of structure functions to Euclidean space offers many 
simplifications. First of all, it gives an access directly to the subject of pQCD, the Q 2 - 
evolution, leaving only one undefined parameter to the measurement. This parameter is 
the absolute normalization of the moment at some arbitrary scale. This scale can be chosen 
in the kinematic region where approximations involved in calculations are fulfilled. The 
second, the calculation of higher order corrections by means of effective series summation 
techniques are less involved in the moment space Q. Obviously, there are also difficul- 
ties: nucleon structure functions cannot be measured down to x = 0. Therefore, lowest 
moments always carry a dose of uncertainty. However, higher moments are well defined 
experimentally and they are subject to precise measurements |2[ pi. 

Scaling down in Q 2 we explore a new kinematic region, not yet understood in terms 
of pQCD. It is the region of bound partons, where the interaction with the probe photon 
likely involves more than one single parton. Correlated partons can either interact to col- 
lectively produce hadronic final state or form an excited nucleon state. These multiparton 
correlations are responsible for the confinement and they are the subject of intensive stud- 
ies j|, H ||, 0, |§|, ||, [u| . Main goal of these investigations is to estimate the contribution 
of multiparton correlations (higher twists) and their Q 2 evolution based on experimental 
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data. The theoretical basis of these studies j|, 11 1 is still poorly developed. This creates a 
large amount of phenomenological approaches aimed to catch main features of the problem 
as e.g. the parton-hadron duality phenomena [l^, ^ |5[ |(| 1C, 13]. 



In this article we discuss applicability of the Operator Product Expansion (OPE) to 
poles of the forward Compton scattering amplitude. The separation of the pole contribution 
in the moments allows to put a lower limit on the higher twist contribution in the measured 
moments. Moreover subtraction of the pole contribution might help to the application of 
the constituent quark model picture as proposed in Ref. [14|. 

2. Forward Compton scattering amplitude 

The forward Compton scattering amplitude 7V\ of the virtual photon with four-momentum 
q on the nucleon N is given by |15| |: 

T,a = ^/ d ^ < N,a\T(J^)Jx(0))\N,a > , (2.1) 

where the sum is running over polarization degrees of freedom a, are hadronic currents, 
£ is the light-cone separation between currents and T() stays for the time-ordered product. 
In case of unpolarized scattering off the nucleon the amplitude can be expressed in terms 
of two invariant amplitudes T\ and T 2 : 

^ = -r a + q -$)+w^ + 9^) ( Pa + 9^) ' (2 - 2) 

where q^ = (u, q) (Q 2 = —q 2 ) and are the virtual photon and proton four-momenta, 
respectively. 



Hadronic currents from Eq. 2.1 for particular class of processes realizing through the 
production of an intermediate state particle can be expressed in terms of invariant elastic 
and transition form-factors. The nucleon intermediate state corresponding to the elastic 
current according to Ref. [|l(| can be written as: 

< N',P\j£(0)\N,a >= u p (P ! ){^F 1 (Q 2 ) + i^F 2 {Q 2 )}u a {P) , (2.3) 

where -Fi and F 2 are known Dirac and Pauli form-factors and M is the proton mass. The 
intermediate particle can be also a nucleon resonance. In this case the hadronic current 
can be expressed as following flT^] : 

2 

e 



< Nj, P\eZJ&(p)\Ni, a >= - ^ M (M 2 R - M 2 ) A m , (2.4) 



where j is the resonance spin, Mr is the resonance mass, A m is the resonance helicity 
amplitude, a and (5 are helicity indices, m indicates the virtual photon helicity state with 



respect to the proton helicity defined as in Ref. [16|, and the virtual photon polarization 
four- vectors are given by: 

e± = -^(0,±l,-i,0), e° =1(131,0,0,1/) . (2.5) 
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These two contributions of the amplitude are particular because represent poles on the 
Riemann ^-surface. The rest of the amplitude is an analytic function of v describing the 
continuum of intermediate states. Therefore, the total amplitude can be described as a 
sum of elastic (T E ), resonance (T R ) and continuum (T c ): 

T(u, Q 2 ) = T E (v, Q 2 ) + T R {u, Q 2 ) + T c {u, Q 2 ) . (2.6) 

If we plug the complex amplitude from Eq. |2.6| into the unitarity relation SS^ = 1 (with 
5 = 1 + 2iT) we obtain 

ImT E (v,Q 2 )+ImT R (v,Q 2 ) + ImT c Xv,Q 2 ) = \T E \ 2 + \T R f + \T C \ 2 + 2Re[T R T C] ] , (2.7) 

collecting the same singularities in r.h.s. and l.h.s. and neglecting the interference effects 
we obtain the optical theorem for each individual term separately. 

In this article we are interested in the case of T 2 . We assume that elastic and resonance 
amplitudes are simple poles and can be written in the factorized form: 

uT E,R = R E,R {Q 2 )p E,R {u) ^ (2 g) 

where the residue of the elastic pole is given by 

*•(«•) = <3W>+*frW> , (2 . 9) 

1 + r 

where Sachs form- factors Ge and Gm are related to F\ and F 2 as following: 

G E = F 1 - ^F 2 , G m = F 1 +F 2 . (2.10) 
While resonance residues can be written as 

* H <« 2 > = i%^K>l 2 + l^ + 2^| Sl/2 |=l , (2.11) 



here qR, is the virtual photon three-momentum taken at the resonance pole. 
3. Pole structure of the Compton amplitude 

The Compton amplitude defined in the first quadrant of the complex z/-plane can be an- 
alytically continued in other quadrants by means of the crossing symmetry. The analytic 
continuation can be performed through the following relation between crossing channels: 
Tzi—v) = T2(u) (time reversal) and T 2 (u*) = T 2 *(z^) (77 — ► pp) . The Compton amplitude 
has a number of simple poles and branch cuts shown in Fig. H. All these poles and cuts 
should obey crossing symmetry relations mentioned above. The elastic scattering ampli- 
tude has two poles situated on the real axis at u = ±ve = ±Q 2 /2M and it has the following 
simple structure: 

P E {v) = + . (3.1) 

V E -V V E + V 



1 in the following two sections we will drop Lorenz indices for simplicity 
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Nucleoli resonances also generate poles in 
the Compton amplitude. But these are located 
on unphysical Riemann sheets. In order to sat- 
isfy symmetry relations each resonance should 
have at least four poles, two on each lower and 
upper unphysical Riemann sheets 2 . Neverthe- 
less, we limit ourselves to simple poles, which 
do not carry information about Riemann sheet 
of location 3 : 
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Figure 1: The integration contour C con- 
tains usual Re v branch cuts from ±Q 2 /2M 
up to ±oo and poles of the nucleon excited 
states. Poles are situated on the unphysical 
Riemann sheets (poles indicated with dot are 
located at lower sheet and those with crosses 
are on the upper sheet) and do not appear 
within the contour C. 



V R + V 

(3.3) 



where e is an infinitesimal positive constant in- 
dicating Riemann sheet of influence and the res- 
onance is located at: 



vr 



Q 2 + (M 2 -M 2 ) T r M r 



(3.4) 



2M 2M 
where Tr is its width. One can check that 
Eq. [T^ yields the familiar Breit-Wigner param- 
eterization. 

Unfortunately the simple parameterization from Eq. |3.3| does not respect the second 
crossing symmetry property: T^i^*) = l^iy) and the amplitude becomes imaginary on the 
real axis even in the region \u\ < Q 2 /2M where no branch cut is present. We can adjust 
this by brute force method multiplying the +ie term in Eq. by: 

|-1,H>QV2M 

[ 1, \v\ < Q 2 /2M . ( 1 

This picture is in contrast with Refs. (R| where only poles on the lower unphysical 
sheet are predicted. However, we find that the absence of poles on the upper unphysical 
sheet would generate a disbalance between lower and upper semi-planes of the Compton 
amplitude on the physical sheet. This disbalance would result in an imaginary part of the 
discontinuity across the branch cut and therefore would violate the optical theorem: 



Disc{T R ) 
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T R {y + ie) £ C ^ 2txW r {v) £ R , (3.6) 



where W R is the resonance part of hadronic tensor. 
4. Low Energy Expansion 

The standard OPE derivation for DIS proceeds through an expansion of the dispersion 

integral in the series of v in unphysic al kinematic domain \v\ < Q 2 /2M, which we call here 
2 actually poles are present on many lower and many upper unphysical sheets generated by different 
decay channels 

3 in principle the pole structure of the resonance amplitude is more complex: poles should "remember" 
the Riemann sheet which they belong to. This can be parameterized in various forms. One of the simplest 
solution is the following: 

P» = -i , (3.2) 

\JVB ± V — \JVB — VR 

plus two conjugate poles at v R , here ± signs refer to the principal and secondary square root values 
indicating therefore correct Riemann sheet. Howeyeg, in this study we are interested in the amplitude at 
small v « ve where the difference between two cases is not so important but the simplification of using 




Figure 2: The imaginary part of the resonance shape distribution: the solid line shows the shape 
of Eq. |3.3| (Breit-Wigner) and the dashed line represents Eq. 3^2. Note that the dotted line goes to 
zero at the threshold. 



"low energy" regime. The Cauchy integral over the contour C shown in Fig. p] is given by: 

12 ,J\J,,I 



i ,T^yw (41) 

2m J c v ' — v 



Assuming that the integral over the circle of radius R — ► oo vanishes and taking dis- 
continuities across left and right branch cuts, using also the time reversal invariance 
T(Q 2 , -v) = T(Q 2 , v) we obtain: 

T(Q 2 ,u')du' _ [°° Disc(T(Q 2 , v'))dv' [~ UB Disc(T{Q 2 ,v'))dv' 



c V - v 



2i I I — !— + —L_ \Disc(T(Q 2 , v'))dv' . (4.2) 



I'E 



v' — v v' + V 

From the other hand the optical theorem states that: 

Disc(T(Q 2 , v)) = Im{T(Q 2 , v + is)) = 2ttW(Q 2 , v) . (4.3) 

Therefore, substituting the discontinuity DisciT (Q 2 , u)) by Eq. 4.2 we can rewrite the 
dispersion relation in the following form: 

T(QV) =4 r -^- 2 {w{Q 2 ,A . (4.4) 
Choosing \v\ < Q 2 /2M one expands this into the geometrical series: 

oo , 

r 1 ~i ti 

T(QV) = 4 Y. M n {Q 2 ) - , (4.5) 

n=0,even 

where 

M n (Q 2 )= [ dx'x' n+1 W(Q 2 ,x') . (4.6) 
Jo 



- 5 - 



In turn, the OPE of the product of two hadronic currents separated by a small light- 
cone distance is given by 



T(Q 2 ,u)=^Crn(Q 2 ^ 2 )Ol(^)(l) (± 



t/2-1 



(4.7) 



and making the standard correspondence between terms with equal powers of 1 jx we obtain 
usual DIS twist expansion: 



T=2 



(4.8) 



where expansion coefficients C Tjn can be analytically calculated in pQCD and local opera- 
tors 01 can be obtained from Lattice simulations. 



In order to apply this expression in practice we need to rewrite Eq. 4.5 for a measured 
structure function. For instance, it can be done for moments of the proton structure 
function F<i = vW 2 extracted from the data in Ref. n-th moment of the structure 
function F 2 is defined as: 



M n (Q 2 ) = f dxx n F 2 (x,Q 2 ) , 
J 



and therefore one obtains the following expansion: 



2 \ Ur-2) 



(4.9) 



(4.10) 



Mn(Q 2 ) = Enr(V,Q 2 )O nT (n)( ^ 

r=2k VV ' 

where k = 1,2,. ..,00, n = 2,4...,oo 4 , fi is a reference scale, O nT (/x) is the reduced matrix 
element of the local operators with definite spin n and twist r (dimension minus spin), 
related to the non-perturbative structure of the target. E nT (fi,Q 2 ) is a dimensionless 
coefficient function describing the small distance behavior, which can be perturbatively 
expressed as a power expansion of the running coupling constant a s (Q 2 ). 

However, we noticed that also the pole (elastic and resonance) amplitude contributing 
to the l.h.s of Eq. 4.13 can be expanded in the power series of v in the region \ v\ < Q 2 /2M: 

00 ~ 00 -. 

vTf(v,Q 2 ) = 2R E (Q 2 ) £ pq" = V(Q 2 ) £ [±1" (4.11) 



n=l,odd 



n=0,even 



and 



00 

uT 2 R (u,Q 2 ) = 4R R (Q 2 ) £ 



n=l,oda! 



4 \ r 

cos (n^fl) = -Izfll 2^ I" 

n=0,etien 



XR\ 
X 



cos ((n + l)0i2) , 
(4.12) 



where B = atan Q2+ ^^ M2) is the phase of z/ fl . 

4 n = relation is meaningless because Regge theory predicts that the structure function behaves as 1/x 
at x — ► and therefore n = integral in Eq. 4.9 would be diverging. 
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In case of T% amplitude the moment expansion will give 



uT 2 (Q 2 ,u) = - V M n (Q 2 ) 



n=0,even 



L X 



(4.13) 



We can rewrite Eq. 4.13 as following: 



vT°{v,Q 2 ) = - V M n (Q 2 
x ' 



1 



,r 



vTfiy^-vT^Q 



(4.14) 



n=0,etien 

therefore the continuum part of the Compton amplitude is given by 

A 00 1 1 

*/T|>,Q 2 ) = - ^ y n {Af n (Q 2 )--^(Q 2 )- J R«(Q 2 )|x i? r +1 cos((n + l)^)} 



n=0,ei>en 



or in terms of measured moments 5 it can be rewritten 



(4.15) 



1 



M n c (Q 2 ) = M n (Q 2 ) - -Ff(Q 2 ) - R R {Q 2 )\x R \ n - 1 cos ((n - 1)^) , (4.16) 

where the last two terms can be estimated using phenomenological data on elastic and 
resonance transition form- factors. 



5. Results 

First of all one can notice that the elastic contribution comes out to be a factor of two 
smaller than in all previous works (see for example Ref. 0). Previously, the dispersion 
relation method was applied uniformly to both the elastic and inelastic channels. This 
led to a confusion in two aspects: 1) the discontinuity across a pole is not a well defined 
object from the mathematical point of view; 2) the value of the integral along the semi- 
circle around the pole is not vanishing. One can prove this statement by evaluating the 
dispersion relation for the elastic amplitude alone. If the elastic amplitude is given by: 



T(u) 



1 



1 



+ 



(5.1) 



then applying the dispersion relation from Eq. 4.1 using the discontinuity across the pole 
one finds: 



1 



■ + ■ 



1 



h A, 

VE — V VE + V 7T Jo 



1 



1 



V — V 



x lmT(v')dv' 



1 



1 



VE-V VE + V 



(5.2) 



Therefore, the evaluation of the discontinuity across the elastic pole in the dispertion 
relation leads to the amplitude which is twice larger than the original one. The comparison 
of moments obtained in the present article and those from Ref. [jTfl is shown in the Fig. || 
together with various pQCD calculations. 



5 notice that the definition in Eq. 4.£ is different from that of the Ref. El and we replace here n with 
n — 2 for consistency 
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Figure 3: Total Nachtmann moments of the proton structure function F2 from Ref. M obtained 



by adding the elastic contribution according to Eq. 4.16 (solid triangles) and "standard" formula 
(empty circles). The lines represent pQCD calculation fitted to the data at largest Q 2 : solid line - 
LO, dashed line - NLO, dotted line - NNLO, dot-dhashed - NLO+SGR from Ref. @. 



Using the relation from Eq. 4.1£ we estimated the resonance contribution to the proton 
structure function moments. For the sake of simplicity we neglected the longitudinal cou- 
plings S1/2 of nucleon resonances. For this exploratory work this is a good approximation, 
provided that those couplings are typically small. We used phenomenological parameteri- 



zations of the nucleon resonance helicity amplitudes from Ref. and Ref. [20]. We have 
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Figure 4: Inelastic Nachtmann moments of the proton structure function F% from Ref. B: solid 
line represents the resonance contribution with form- factors of Ref. [pl|; dashed line shows the 
resonance contribution using form- factors from Ref. p3 . 



taken into account 9 (20) well established resonances from the parameterization J2(| (Jl9|]) 



whose masses and widths were obtained from Ref. [21]. In order to compare this to the 



data extracted in Ref. || we had to convert Conwell-Norton moments into Nachtmann 
ones. In the region M 2 /Q 2 < 1 this can be done by expanding Nachtmann moments in a 
power series of M 2 /Q 2 as in the Ref. |^ ]. Results of the comparison are shown in Figs. [| 
and Fig. ||. 
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Figure 5: The total higher twist contribution and the elastic+resonance contribution to the 
F% structure function moments relative to its leading twist part taken from Ref. ||: solid line - 
the total higher twist contribution; dashed line - the resonance contribution using resonance form- 
factors from Ref. |20) ; dotted line - the elastic+resonance contribution using resonance form-factors 
from Ref. [f(i|. 

6. Conclusions 

We calculated the contribution of the nucleon elastic pole and its excited states into mo- 
ments of the nucleon structure function F%. Due to formation of an intermediate bound 
state this contribution cannot be described by an incoherent elastic scattering of the vir- 
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tual photon off a quark (either constituent or current). In terms of OPE this contribution 
belongs to higher twists. This is confirmed by the observed Q 2 -dependence of estimated 
the resonance contribution shown in the Fig. |1 Vice versa, assuming that conceptually 
the resonance contribution as well as the elastic peak contribution should belong to higher 
twists one obtains a stringent test on the large-Q 2 behavior of resonance form-factors. 

We found that the elastic peak contribution is actually twice smaller than it was known 
before. The reduction of the elastic contribution in moments leads to an improvement of 
the duality effect in the low Q 2 region. As you can see in Fig. || the second moment 
becomes almost constant while higher moments are much more flat after the correction 
and leave no room for a significant higher twist contribution (see also Fig. ||). This leads 
to the conclusion that the extraction of nucleon elastic form-factors from DIS pioneered in 
Ref. |?J is not hopeless and the large discrepancies found more recently in Ref. ||| are due 
to missing 0.5 factor. 

In Refs. [^, |H we have seen that different higher twist terms cancel generating the 
well known duality phenomenon. The resonance contribution estimated in this article 
yields coherent positive higher twist in the moments. Additional positive contribution to 
higher twists is given by the elastic scattering. While the negative higher twist which 
largely cancels the resonance and elastic contributions is related to number of exclusive 
channel thresholds (e.g. single pion production threshold). If all elastic, resonance and 
the threshold contributions are subtracted from the moments the remaining part should be 
related to the incoherent scattering off proton constituents regardless Q 2 value. Therefore, 
an estimate of the threshold contribution to the structure function moments would allow 
to study Q 2 -behavior of form-factors of proton constituents. 

Notice that both elastic and resonance higher twist contributions are related to nucleon 
and resonance bound state wave functions. Therefore these terms cannot be described by 
asymptotic freedom methods like Borel summation of the leading twist as series which also 
generates 1/Q 2 corrections. Given large magnitude of elastic and resonance contributions 
in the total moments (in particular for n > 2) the good overlap between renormalon-based 
fits [^4|] and the data should be taken as accidental or related to the duality phenomenon. 

It is clear from Fig. [5| that the resonance contribution itself does not follow DIS scal- 
ing behavior. Therefore, the observed approximate scaling of integrated resonance peaks 
(known as local duality) is due to compensation between rising resonances and decreasing 
background at low Q 2 . 

Large physics program focused on extraction of nucleon resonance form-factors from 
electro- and photo-production data undergoing in Jefferson Lab will allow for precise eval- 
uation of the main higher twist contribution to the nucleon structure functions. 
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